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Critical Slowing Down of Triangular Lattice Spin-3/2 Heisenberg 
Antiferromagnet Li 7 Ru06 via 7 Li NMR 
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Kanako NAKAJIMA 2 , and Hirohiko SATO 2 
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We report 7 Li NMR studies of single crystals of triangular-lattice Heisenberg antiferromagnet 
Li7Ru06- Slow critical divergence with a wide critical region of |T/Tn — 1 1 < T was observed 
in 7 Li nuclear spin-lattice relaxation rate. The slowing down of staggered spin fluctuations was 
analyzed in a renormalized classical region of a two-dimensional triangular-lattice non-linear 
sigma model. A spin stiffness constant was found to reduce to about 20 % from the value in 
a spin-wave approximation. The effect of spin frustration, e.g., Z2 vortex excitations on the 
critical phenomena is suggested. 
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Spin frustration effect of a triangular lattice Heisen- 
berg antiferromagnet has been one of the central issues 
in physics and chemistry of magnetic insulators. The 
ground state is classically a long range ordering state 
with the 120° spin structure, but the quantum mechan- 
ical ground state might be a gapped or gapless quan- 
tum spin liquid. 1-3 Elementary excitations of Z2 vortices, 
topologically stable point defects, are inherent even in 
classical triangular spin systems. 4-7 At finite tempera- 
tures, topological phase transition between the param- 
agnetic states is theoretically predicted from numerical 
simulations. 

Finite-temperature magnetic long range ordering of 
quasi low dimensional antiferromagnets may be ascribed 
to a three dimensional interaction and the anisotropy. 
The low dimensional characteristics are (i) the sup- 
pressed magnetic ordering temperature Tn < T^f F in 
the mean field approximation, (ii) the suppressed Curie- 
Weiss law, i.e., the maximum of uniform spin susceptibil- 
ity at T max due to short range ordering, and (iii) a wide 
critical region in two dimensions. 3 

Two dimensional square lattice Heisenberg antiferro- 
magnet has been intensively studied through the stud- 
ies of high-T c cuprates. 2,3 Our understanding a short 
range ordering has made rapid progress. In the renor- 
malized classical region, an antiferromagnetic correlation 
length and a staggered spin susceptibility diverge expo- 
nentially toward Tn = K, 8 ' 9 which were observed in 
La2Cu04. 10, 11 Such a wide critical region toward Tn = 
K is the notable feature of two dimensional square 
lattice Heisenberg antiferromagnets. Non-linear sigma 
model and field theoretical treatments turned out to be 
the relevant model and to give us powerful methods to 
describe the low energy excitations. They were also ap- 
plied to the non-collinear frustrated magnetic systems 
like the triangular lattice. 12-17 To our knowledge, how- 
ever, the theoretical critical behavior has not been widely 
tested to the actual critical phenomena of triangular spin 
systems. 
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A delafossite-type I^RuOg has been at first reported 
as hexagonal "Lis-^RuOg" . 18 However, a detailed anal- 
ysis on powder 19 and on single crystal 20 revealed that 
the actual composition is "LiyRuOg" and the struc- 
ture is slightly distorted from an ideal hexagonal lattice. 
LiyRuOe is triclinic and has the superlattice structure of 
the Li deficiency in the double Li layers. 20 Figures 1(a) 
and 1(b) show the crystal structure and the triangular 
lattice, respectively. Single Li2RuO@ layer and double Li 
layers are stacked with each other. The Li2RuOg layer 
consists of a triangular lattice of RuOg octahedrons and 
a Li honeycomb lattice. The Li ions occupy 12 sites in 
the RuOg layer and 30 sites in the Li deficient layers in 
unit cell, which are crystallographically inequivalent. For 
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Fig. 1. (a) Schematics of crystal structure of Li7Ru06 and (b) the 
top view of a triangular lattice Ru plane, (c) Inverse magnetic 
susceptibility x _1 an d (d) the magnetic susceptibility x of a 
single crystal Li7RuOg. The solid lines are the best fits by an 
inverse Curie- Weiss law. The arrow indicates Tn ~ 6.5 K. 
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a hexagonal model of LigRuOg, the staggered magnetic 
field from Ru would be cancelled out at the Ru-plane 
Li site. However, for the actual LiyRuOg, since the Li 
ions occupy non-ideal positions deviated from the center 
of the Ru triangle, then the 120° staggered spin fluctu- 
ations can be probed through the Li NMR. Li7RuC>6 is 
a deformed triangular lattice system. The detail of the 
crystal growth and the structure analysis will be pub- 
lished in a separated paper. 20 

Uniform magnetic susceptibility Xa {ot = | and _L 
denote H | plane and H _L plane) shows a slightly 
anisotropic Curie- Weiss behavior \ = C/(T — 6) with 
C = 1.89 emu/mol-Ru and 6 = — 73 K at high temper- 
atures. Figures 1(c) and 1(d) show the inverse magnetic 
susceptibility x _1 up to 300 K and the magnetic suscep- 
tibility x on an enlarged scale for a single crystal, respec- 
tively. Upon cooling, Xa makes a maximum at about 12.5 
K and drops sharply at the Neel temperature Tn ~ 6.5 

19,20 rp-^g k roac l maximum behavior indicates a low di- 
mensional magnet. The Curie constant C is close to the 
value for S = 3/2 and g = 2. For a Heisenberg spin Hamil- 
tonian Y>J nn Si ■ Sj with the z nearest-neighbor exchange 
interaction J nn between Ru ions, the Weiss temperature 
is given by 



(a) 



6 = 



S(S + l) 



(i) 



in the mean field approximation. Using S = 3/2 (Ru 5+ ) 
and z = 6, we estimated the superexchange interaction 
J nn = - 9.7 K. Since T N w 6.5 K and 6 = -73 K, we ob- 
tain the ratio of Tn/|B| « 0.089, which is nearly the same 
as 0.082 for VC1 2 with S = 3/2 (V 4+ ). 21 These are the 
low dimensional characteristics of (i) and (ii). Thus, the 
layered compound Li7RuC>6 is a quasi-two dimensional 
antiferromagnet . 

In this Letter, we report 7 Li NMR studies of single 
crystals of triangular-lattice Heisenberg antiferromag- 
net Li 7 Ru0 6 . We found slow critical divergence with 
a wide critical region of |T/T N — 1| < 7 in 7 Li nuclear 
spin-lattice relaxation rate, which is due to the slowing 
down of staggered spin fluctuations. From the analysis 
by renormalized classical fluctuations, we found the re- 
duction of a spin stiffness constant, possibely due to the 
spin frustration effect. 

Single crystals of LiyRuOe were grown from the mix- 
ture of RuC>2 and a large amount of LiC03 flux heated in 
an oxygen atmosphere at 950 °C. Typically 0.5x0.5x0.01 
mm sized and plate-like single crystals were obtained. 20 
The plane and the vertical axes are confirmed to be the 
ab plane and the c axis of the quasi-hexagonal lattice, 
respectively. X-ray diffraction patterns for the powdered 
samples indicated the samples in a single phase. We per- 
formed 7 Li (nuclear spin 7 = 3/2 and nuclear gyromag- 
netic ratio j n /2iT = 16.546 MHz/T) NMR spin-echo mea- 
surements at H = 7.48414 T for the samples in which 
the several single crystal pieces were put on a plate. 
The NMR frequency spectra were obtained from Fourier- 
transformed spin-echoes or summing them at several fre- 
quencies and frequency-swept spin-echo intensity. The 
7 Li nuclear spin-lattice relaxation curves were measured 
by an inversion recovery technique and the relaxation 
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Fig. 2. (a) 7 Li NMR frequency spectra of the single crystals of 
Li7Ru06 at H _L plane (solid curves and closed circles) and 
at H || plane (open circles). The vertical line at 123.8343 MHz 
indicates the 7 Li NMR spectrum peak of LiCktg for a reference 
of zero shift. The inset shows the best fits by three Gaussian 
functions, (b) Knight shift 7 K plotted against the bulk magnetic 
susceptibility \ a with temperature as an implicit parameter. 



times Ti 's were obtained from the fits by a stretched ex- 
ponential function of exp{-(t/Ti)^} [t is the time after 
the inversion pulse). 

Figure 2(a) shows 7 Li NMR frequency spectra of the 
single crystals of Li 7 Ru0 6 at H _L plane (solid curves 
and closed circles) and at H \\ plane (open circles). At 
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Fig. 3. (a) Temperature dependences of 7 Li nuclear spin-lattice 
relaxation rates 1/Ti of SI (closed circles) and S2 (open trian- 
gles) at H _L plane and T > T^. The inset figure shows semi-log 
plots of the temperature dependences of the stretched exponents 
(5. (b) Log-log plots of 7 Li nuclear spin-lattice relaxation rates 
1/Ti at H _L plane against the reduced temperature (T— Tn)/Tn. 
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T > 200 K, a single sharp 7 Li NMR spectrum and no 
quadrupole splits were observed. Upon cooling, the NMR 
spectrum shifts to higher frequency side and a weak sig- 
nal separates and shifts more largely. At 4.2 K, both 
NMR spectra at H _L plane and H || plane are broad- 
ened nearly symmetrically. The featureless symmetric 
broadening below Tn indicates the emergence of internal 
magnetic field of an incommensurate staggered moments 
along the ab plane and the c axis. 

The inset shows three Gaussian functions fit to the 
NMR spectrum at H _L plane and T =19.8 K. A sharp 
strong peak and a broad higher frequency peak are de- 
noted by SI and S2, respectively. S3 just adjusts the foot 
of the spectrum, whose Knight shift is ~T00 ppm and 
nearly independent of temperature. From the NMR in- 
tensity, the strong peak SI and weak S2 can be assigned 
to the Li sites in the double Li layers and in the RuOe tri- 
angle lattice layer, respectively. The assignment at H \\ 
plane at T > T N not shown in Fig. 2(a) was less clear, 
since at least 4 peaks were observed. The lowest and 
highest frequency peaks are assigned to SI and S2, rc- 
sepctively. 

Figure 2(b) shows 7 Li Knight shifts of SI and S2 plot- 
ted against the bulk magnetic susceptibility \a with tem- 
perature as an implicit parameter. The bulk magnetic 
susceptibility Xa is expressed by the sum of a tempera- 
ture dependent spin susceptibility x sp in, the Van Vleck 
orbital susceptibility xvv, and the diamagnetic suscep- 
tibility of inner core electrons Xdia- The 7 Li Knight shift 
is expressed by K a = A a (i)x sp in,a/N A HB+K dia (i =1 
and 2 denote SI and S2, Na is Avogadro's number, and 
/iB is the Bohr magneton). From the K — \ plots in 
Fig. 2, we obtained positive hyperfine coupling constants 
of A,|(l) w A ± (l) = 0.45 kOe/^B, A\\(2) = 2.0 kOe/^ B 
and Aj_(2) = 1.2 kOc/^ B . The positive A(l) and A(2) 
indicate the predominant role of a transferred hyperfine 
coupling. 

Figure 3(a) shows temperature dependences of 7 Li nu- 
clear spin-lattice relaxation rates 1/Ti of SI and S2 at 
H _L plane and T > Tn. The inset figure shows semi- log 
plots of the temperature dependences of the stretched ex- 
ponents j3. At high temperatures, 1/Ti levels off, which 
indicates the exchange narrowing limit. At low temper- 
atures, 1/Ti shows the divergence behavior due to the 
critical slowing down of the staggered spin fluctuations. 

In general, the high temperature limit of 1/Ti m the 
exchange narrowing is expressed by 22 
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The hyperfine field at the in-plane Li results from the 3 
nearest neighbor Ru spins. The number of the nearest 
neighbor exchange coupled Ru spins is 6. Putting z — 
6, S = 3/2, A||(2) = 2.0 kOc/^ B and J nn = -9.7 K for 
eqs. (2) and (3), we obtain 1/Tioo = 340 s _1 , which is 
the same order of magnitude but slightly smaller than the 
actual 1/Ti ~ 400 and 600 s^ 1 for SI and S2 above 200 
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Fig. 4. Log-log plots of the nuclear spin-lattice relaxation rate 
(l/Ti)/(l/Tioo) or the continuous-wave (cw) NMR lincwidth 
A/Aoo against the reduced temperature (T — Tn)/Tn for the 
critical phenomena of CuO, 24 FeF 2 . 25 K 2 MnF 4 , 26 K 2 NiF 4 , 27 
La2Cu04, 10 and Li7Ru06 (SI). The circles, squares and trian- 
gles are three dimensional magnets, quasi two dimensional square 
lattices, and a triangular lattice, respectively. The solid and open 
symbols are the data by spin-echo and continuous-wave (cw) 
measurements, respectively. 



K, respectively. The experimental value of Ti depends 
on the fitting function more or less. Using the fitting 
function of a double exponential function, we obtained a 
long component of 1/Ti ~ 277 and 435 s _1 for SI and 
S2, respectively, being the same order of magnitude of 
the estimated 1/Tioo. 

Figure 3(b) shows log-log plots of 7 (1/Ti) against the 
reduced temperature (T — Tn)/Tn. The critical diver- 
gence of 1/Ti starts from |T/T N - 1| ~ 7. Although such 
a wide critical region is not usually regarded as a three 
dimensional critical region, we tried to apply a power law 
of |T N /(T - T N )|" and then obtained n = 0.45. This is 
close to the mean field value. 23 

Figure 4 shows log-log plots of the nuclear spin-lattice 
relaxation rate (1/Ti)/(1/Ti 00 ) or the continuous-wave 
(cw) NMR linewidth A/A^ against the reduced tem- 
perature e =(T — Tn)/Tn for the critical phenomena 
of CuO, 24 FeF 2 , 25 K 2 MnF 4 , 26 K 2 NiF 4 , 27 La 2 Cu0 4 , 10 
and Li7RuC>6. The onset temperatures of the increase in 
Tioo/Ti and linewidth A/A^ are the beginning of the 
critical slowing down of relevant spin fluctuations and 
indicate the width of the individual critical region. For 
a narrow critical region of |e| < 0.1 of CuO and FeF 2 , 
the critical exponent of three dimensional critical slow- 
ing down was estimated. 24, 25 The wide critical region of 
|e| < 1.0 of La 2 Cu0 4 was successfully understood by a 
two dimensional short range ordering effect in the renor- 
malized classical region of the square lattice. 10 The wide 
critical regions of square lattices K 2 MnF 4 and K 2 NiF 4 
were also understood by the renormalized classical fluc- 
tuations. 28, 29 The critical region of |e| < 7 of LiyRuOg 
is wider than those of the two dimensional square lat- 
tice spin systems. One should note the another striking 
feature of the slow divergence of 1/Ti of Li7RuOg. The 
wide but slow critical divergence of 1/Ti characterizes 
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Fig. 5. Semi-logarithmic 1/TiT 3 against 10/T for Li 7 Ru0 6 (SI) 
at If 1 plane. The straight line is the fitting function of eq. (5). 



the critical phenomenon of L^RuOe- 

For the frustrated quantum antiferromagnets, the non- 
linear sigma model description tells us the magnetic cor- 
relation length 12 - 14 ' 15 



£ oc -^=cxp(47rp Q /T) 

with a spin stiffness constant p a (a = _L and | 
and the nuclear spin-lattice relaxation rate 14 ' 15 

±- cx r 3 cxp(47rp Q /r). 
J-l 

The spin stiffness constant p± is given by 

1.51|J|, 



P± = ^Zj.SC\J\ 



(4) 
plane) 

(5) 
(6) 



where a renormalization factor Z±_ is estimated by a spin- 
wave approximation and 1/5 expansion. 16, 17 Equations 
(4) and (5) are applicable to the low temperature states 
at T < 2tt ( o q . 14 ' 15 

The NMR linewidth can serve as a probe of the mag- 
netic correlation length £ in real space. The linewidth of 
the S2 NMR spectrum in Fig. 2(a) increased rapidly at 70 
K upon cooling, which did not scale with the Knight shift 
K but was similar to the divergence in \/T\. Thus, we 
may assume the temperature dependent £ as in eq. (4). 

Figure 5 shows semi-logarithmic \ jT{T z against 10/T 
for LI7RUO6 at H _L plane. The straight line is the fit- 
ting result of eq. (5). From the slope of the line at lower 
temperatures T < 15 K, we estimated the exchange in- 
teraction |J| = 2.1 K. Since 2np± = 20 K, the fit range 
is justified a posteriori. The interaction \J nn \ = 9.7 K 
from the Weiss temperature reduces to |J| = 2.1 K 
at low temperatures. Since J nn should not strongly de- 
pend on temperature in the range of 4.2 K < T < 300 K, 
the alternative reduction should trace back to the renor- 
malization factor Zj_ and the stiffness constant p±. Z±_ 
should reduce to about 20 % from the value in the spin- 
wave approximation. This reduction might be the effect 
of the spin frustration, e.g., the topological Z 2 vortex ex- 
citations on the spin correlation, which agrees with the 
numerical simulations. 7 The reduced Z± and p± are just 
to put the slow divergence of 1 jT\ in another way. 

The wide critical region has been observed for the 
other 5 = 3/2 triangular lattice systems such as VCI2, 30 
HCr0 2 , 31 and LiCr0 2 . 31 ' 32 It is, however, less clear 



whether the slow critical divergence is ubiquitous in the 
triangular lattices. 

In conclusion, we found slow critical divergence with 
a wide critical region of |T/T N - 1| < 7 in 7 Li nuclear 
spin-lattice relaxation rate for the triangular-lattice anti- 
ferromagnet L^RuOe . We applied renormalized classical 
staggered spin fluctuations to the slow divergence and 
then obtained the reduction of a spin stiffness constant, 
suggesting the spin frustration effect. 
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